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Abstract 

In the gauge-mediated supersymmetry breaking models with scalar tau as 
the next-to-lightest supersymmetric particle, a scalar lepton may decay dom- 
inantly into its superpartner, tau lepton, and the lightest scalar tau particle 
through I £r±f{F. We give detailed formulas for the three-body decay am- 
plitudes and the polarization asymmetry of the outgoing tau lepton . We find 
that the tau polarizations are sensitive to the model parameters such as the 
stau mixing angle, the neutralino to slepton mass ratio and the neutralino 
mixing effect. 
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I. INTRODUCTION 



Supersymmetry provides a natural solution to the hierarchy problem that associates with 
weak scale physics and Planck scale physics If nature is supersymmetric, it must spon- 
taneously break supersymmetry and transmit the message to low energy effective models. 
Two approaches are widely considered to satisfy the message transmitting strategy. One is 
the gravity mediation approach which has the supersymmetry breaking scale a/F about of 
order 10^^ GeV [0,^. The gravitino mass is typically much larger than those of gauginos 
and sfermions under this framework. Due to the largeness of the gravitino interaction 
is also too weak to play important roles in collider phenomenology. Another approach uses 
gauge interactions to mediate supersymmetry breaking effect [|,||. In the gauge- mediated 
supersymmetry breaking (GMSB) scenarios, the supersymmetry breaking scale a/F is suf- 
ficiently small so that the gravitino G is always the lightest supersymmetric particle (LSP) 
with mass 
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nrin 



,n - = 2.37 X 10-^ X eV, (1) 

where Mpi denotes the reduced Planck scale ~ 2.4 x 10^^ GeV. Since the soft masses of 
other superparticles are associated with the weak scale in the GMSB models, within typical 
range of \fF, the gravitino is far lighter than other superparticles and becomes the LSP. 

It is well known that different scenarios for the next-to-lightest supersymmetric particle 
(NLSP) in the GMSB models could be crucial in discovering supersymmetry signals in 
colliders and may lead to different phenomenologies. In the neutralino NLSP scenarios, the 
lightest neutralino A^i decays into a photon and a gravitino which carries missing energy and 
escapes the detectors 0. If the SUSY breaking scale \fF is less than a few thousand TeV, 
the neutralino decay A^^i — ^ 76* would be prompt enough to occur within the detector. One 
can therefore use inclusive 77+ -^t + X signals in detecting supersymmetry. In the slepton 
NLSP scenarios, the lightest slepton can promptly decay into its partner and a gravitino 
with a sufficiently small scale \fF 0. Depending on the scale \fF, the slepton NLSP may 
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live long enough to leave tracks and kinks in detectors. In general, other NLSP scenarios are 
also plausible depending upon the various constructions of super symmetry breaking models. 

Scalar leptons can be easily found at future linear e^e~ colliders if it is kinematically 
allowed ||^. By measuring their masses and decay distributions one could determine some 
of the MSSM parameters. The phenomenology of sleptons in future linear colliders then 
deserves detailed study. 

In this paper we concentrate on the GMSB scenarios in which the lighter scalar tau 
lepton (stau) fi is the NLSP. In general, we can write down the stau mass eigenstates fi 
and f2 in terms of the stau mixing angle Of 

fi = cos ^^ff^, + sin 6*^7-/2 (2) 
f2 = - sin 6ffL + cos OfTR, (3) 

where 6f ranges from < < vr. Since the mixing is proportional to fermion Yukawa 
coupling, due to the small mixing efffects for scalar electrons (selectrons) and for scalar 
muons (smuons) , the lightest selectron ci and the lightest smuon /ii are nearly right-handed 
sleptons with almost equal masses. We will use and /i/j to represent the lighter selectron 
and smuon respectively in this paper. 

fi can decay into a tau lepton and a gravitino G in the detectors if the SUSY scale 
a/F is not too large. If the scale a/F is sufficiently large, the stau NLSPs may live long 
enough to leave tracks or kinks in the collider detectors. In many of the GMSB models, 
the lighter selectron and the lighter smuon are merely right-handed sleptons with nearly 
degenerate masses. Unlike fi, they can decay not only through the two-body processes 
cr eG and ftji fiG, but also through the processes such as — > iNi, uiu^fi, 
and Ir — > irfi depending upon the model parameters. For instance, Ir — * irfi will be 
kinematically forbidden if the mass difference [mi^ — m^J is smaller than the tau mass rrir- 
This could happen in the so-called slepton co-NLSP scenarios 0. 

However, in the stau NLSP scenarios, there exists some parameter space where the three- 
body processes Er — > irfi could dominate over the two-body processes Ir —>■ iG and leads 
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to different signatures in colliders |T0[ 



If the decay fi tG talkes place outside the detector, then one can directly identify 
the stau tracks or decay kinks. Thus in the e^e" colliders, the pair-produced sleptons may 
decay in collider detectors through e^e~ — > ^/j^/j — > £"'"£~rrfifi. In the e~e~ colliders, right- 
handed sleptons are produced through the t-channel processes and then decay as e^e^ — >■ 
^ij^ij — i^i^TTTifi. It is also found that e/j and /i/j could have macroscopic decay lengths 
that are smaller than the dimensions of typical detectors if the mass difference l^^i^ — — 
vnf^ I is suitably small. If this is the case, the secondary vertices of the three-body processes 
could be resolved and each fi particle can be identified as coming from one of the pair- 
produced selectrons or smuons. In this case one could measure the polarizations of the r 
leptons. The r polarization provides another observable which could help to pin down the 
stau mixing angle and/or the ratio of the neutrilino and the slepton masses. 

However, there might be the possibility that the two-body and three-body decay modes 
are of comparable amplitudes and can therefore decay in the detector. In this case we might 
have the decay process e^e+ ^r^r ~^ i^i'^r^r^ ^ i^i'^r'^r^ ^ ^, i.e., two tau leptons 
with opposite charges in the final state. One of the tau lepton is coming from the three-body 
decay of the slepton therefore its kinetic energy is less then —171^- — m^J. The other tau 
lepton is from the two-body decay of the scalar tau and hence carries larger momentum and 
is much harder than the tau lepton originated from the three-body decay. In this interesting 
case we can easily identify the softer tau lepton from the three-body decay from ipi- 

Motivated by the above observation, we calculate the polarizations for the tau leptons 
in the three-body decays £r ir^f^ in the paper. It is expected that the tau polarizations 
should depend on the stau mixing angle and the mass ratio m^ymg^. In the section two of 
the paper, we set the notations and give formulas for the polarizations of the tau leptons in 
the three-body slepton decays. Numerical results are presented and discussed in the section 
three. In the section four, we make our conclusion. 
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II. THREE-BODY SLEPTON DECAYS 



As mentioned in the previous section, the three-body decays — > Irfi could dominate 
over the two-body decays 1^ in some parameter space if the SUSY breaking scale \/F 

is sufficiently large. We restrict on the scenario where the two-body processes cr ^ e + Nj, 
are forbidden, i.e., we assume the following mass relation: 



me-\-mf[. > rriej^ > ■me + rrir + nif^ . 



(4) 



Therefore the only competing decay processes for cr and fiR are Ir — >■ ir^fi". 

Because of the invariance of charge conjugation and parity transformation, we can write 
down the following identities for the selectron decays: 



(5) 
(6) 



where L (R) denotes the left (right) helicity for the tau leptons in the decay processes. 
Therefore we will only discuss the decays of e^. 

At tree level, the selectron cr decays through the propagation of bino B and wino W^. 
Since bino B, wino W^, and the two neutral higgsinos hi, h2 can mix and form the neutralino 
mass eigenstates Ni, it is appropriate to do calculation in terms of the neutralino parameters. 
In the basis of -01 = {B, W^, ih\, ih!^), the neutrahno states and the associated diagonal mass 
matrix are given by 



Md = V+MV. 



(7) 
(8) 



Here V denotes the transforming matrix which diagonalizes the neutralino mass matrix M 
( . ...... .\ 



M 



m\ mz sin cos /3 —mz sin Qy/ sin /3 

m2 —mz cos Qw cos /3 mz cos Qw sin /3 

mz sin cos /3 —mz cos Q-^ cos /? — 

—mz sin Q-\j^ sin /3 mz cos sin /3 — 



(9) 



where mz is the mass of Z boson, mi denotes the soft breaking mass of bino m2 is the 
soft breaking mass of wino W^, 9w denotes the weak Weinberg angle and tan/3 is the ratio 
between the vacuum expectation values of the two Higgs doublets in the model, and n is 
the super symmetric Higgs mass parameter. 

In our study, the tau polarizations are measured in the lab frame. When the collider 
beam energy is threshold to produce the slectrons in pairs, the selectrons may travel with 
velocity v before decay and thus experience a boost which is characterized by the boost 



factor 7 = (f/c)/yl — {v/cY. Fig. |l| shows the selectron boost direction in the lab frame. 
The moving direction of the outgoing tau lepton is assigned the Z-direction and the XZ 
plane is spanned by the outgoing electron and the tau lepton. By summing all tree-level 
contributions for the decay c^t^Ti , the transition amplitudes |«ni^^-)p for different 

helicities (L, R) of the tau final states are given by: 

j=i ly '"tvjJ 



+ CijPej^vE + K {sin (p sin sin - + sin0 cos cos ^ ~ "^^^ sin (10) 

Nm^L)!' = I E f.f^'^^ iic^^E^ny'E^ - c;^m^^VETK) sin | 
j=i ly ''t-ffj ^ 

9 6 6 

+ CijPej^y/E — K{sin sin if cos - — sin cos sin - + i cos cos (11) 

where 

g2 = m?^ - 2PEe« - 2PPg^ sin cos(^ - ^) , (12) 

{Eej^,Pejf) and {E,K) denote the (energy, momentum) pairs of the right-handed selectron 
and the outgoing tau lepton in the lab frame respectively, P is the outgoing electron energy, 
rrie^ and m^y denote the masses of the right-handed selectron and the jth neutralino re- 
spectively, Qe represents the electric coupling of electron, and the coefficients Cij,C2j, and 



Vrj are given by ITT 



Z sm o\Y Z cos o\Y 



cos 

Cij = XrVsj cos6^ — V2QeVjij sinOf (15) 
C2j = XrVsj sin Of + V2QeVLj cos Of. (16) 

Here is the coupling strength for the tau Yukawa couphng term XrHiLrRr- Similarly, 

the decay process e~^T~^fi has the transition amplitudes 

4 2Q V* \/E~ 6 
)!' = I E Z^^'^^ ^ {{-Cl^m^^VE^ + C2,E,,VeTK) cos - 
j=i {Q ^f,.) ^ 

Q 

+ C2jPejiVE + /^(sin sin sin - + sin cos cos 2 ^ cos sin 2)}!^ i^'^) 
I'^U" = I E fj^^'^^ i(-Ch^N.^^^+^ + C2,E,,VE^) sin I 

j=l W -^N,> ^ 



+ C2jPeii\/ E — -^( sin sin cos - — sin0 cos sin 2 + ^ cos cos 2)j'l^- (1^) 

As easily seen from Eqs. (p!0| - |l8D , the amplitudes do depend upon the selectron boost angles 
and the selectron momentum P^j^. 

Taking into account the selectron boost effect, the decays rate Tl^r of different tau lepton 
helicities L or R for cr is obtained as 

r(L,ij) = j d(j)d(pF{(j),(p)T(^L,R){'i^^^)^ (19) 



with 



_ r dEd[cose] KP\im^L,R)? ^20) 

J GATT^Eej^ Eej^ — E + K COS 6 — Peji sin (p^siiup sin 9 + COS if COS 9) 
p ^ + ml- ml^ - 2EEe^ + 2Pe^K sin cos 

2[i?ejj — E + Kcos9 — Pgfl sin 0(sin (psin9 + cos (p cos 9)] ' 

where -F(0, ip) is the probability density function for selectron boosting in the (0, (p) direction. 



E is the energy for the outcoming tau lepton, K = ^ [E"^ — ml) denotes the momentum of 
the tau lepton, P is the energy of the outcoming electron, and 9 is the angle between the 
tau and the electron. 

Naively, the decay rates with selectron boost effects may differ to those without selectron 
boost effects to the order of 0(7), i.e. 0{Pej^/Eej^). On the other hand, the spinless nature 
of selectrons leads to the following symmetry for -F(0, (p): 



= n±^)=F{^- <\>, = F(7r - 0, TT ± ^). (22) 

Thus all terms that are linear in cos or sin (yj should be cancelled after integrating out the 
ip dependence. Consequently, by observing that 

^Erest ± Krest = \/E±K{\T ^7 COS siu + 0(7^)} (23) 

where Eyest (Krest) denotes the energy (momentum) of the outgoing tau lepton in the se- 
lectron rest frame, the decay rates in Eq.(^) are actually differ to those without selectron 
boost effects to the order of 0(7^). Therefore, we simply do all the calculations in the 
selectron rest frame in the rest of the paper. 



III. TAU LEPTON POLARIZATION 

The tau lepton polarization asymmetry can be defined as 

P. . <^n>-<^^>, (24) 
^ <Tn> + <TL> ^ ^ 

We can easily calculate the asymmetry by using the formulas presented in the previous 
section. From Eq.s (^ and (^ the CP invariance leads to the following relations for the tau 
polarization: 

P.(e+ ^ e+r+fr) = -^r(e^ - e-r-f+) 

P,(e+ e+r-f+) = -P.(g^ ^ g-r+ff) (25) 

Therefore we calculate only Pr(ej^ — ^ e~^T~^fi ) and PrifiR ~^ e"'"r~f]^) and give the numerical 
results in this section. 

Generally, the transistion amplitudes and the associated decay rates presented in the 
previous section are sensitive to the neutralino masses and mixing effects. For instance, by 
assuming the dominance of A^^i propagating in the tree-level decay processes in the Bino-like 
lightest neutralino scenario, the amplitude functions in the limit A = (mg^ — << mg^j 
are obtained approximately as 
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(26) 
(27) 
(28) 
(29) 



From Eq.s ( p^ - |29|) , it is easily seen that the polarization ratio Pr is determined by the mass 

ratio mj^/nien ■ In the purely right-handed fi limit (sin6'f = 1), Pr becomes 

^ ^!UE(E^--,nrKA-Er ^^^^^ 
St dE eJe^^^(A - B)2 



Eq. (|30|) shows that the Pr ratios always have opposite signs for r+ and r~ production in 
the selectron decays, i.e., the — >• e~^T^fi process tends to have right-handed tau leptons 
in the final states and the e^T^f{ process tends to produce left-handed r"'"'s in this 

limit. 

The tau-lepton polarizations as the functions of the mass difference A with different 
cos6f values are illustrated in Fig. ^ (a) and (b) where the mass ratio = nij^/m^^ = 2 is 
kept fixed. Although it was pointed out in |]TD[ that when the mass differences between fi 
and cr and fiR are less than 10 GeV, the range of cos^f will be within 0.1 < | cos^^f | < 0.3 
, we find that if we assume that mg^ ~ 200 GeV the mass differences between fi and cr 
and fiR can be up to 20 GeV while keeping 0.1 < | cos^f | < 0.3. As shown in the plot (a), 
P(r+) = Pr{(i\i e+r+ff) decreases as A increases. For a given A value, P{t^) always 
has larger values for larger cos^f values. This could be explained by the fact that when 
cos^f increases, the amplitude lini^-jp always increases faster than |iin^-|p does when the 
mixing angle are within the range 0.1 < | cos^^ | < 0.3 . Although not shown in the figure, 
from Eq.s (^,^) one can easily see that when A is kept fixed, a larger cosdf will give 
smaller polarization P{t^) = Pr{e^ —>■ e'^r^fi). 

As shown in Fig. 0(b), the ratio between P{t^) and P{t^) is a slow varying function 
of A with fixed = 2 for most of the cos6'f range except for cos^f = —0.3. This is due 
to the complex interplay between the factors \/E + K cos 9f and \J E — sin Qf. When the 



mass difference A is small, ^/E + K and ^/E — K are comparable and both and 
lilTL^^P are getting smaller as cos6'f changing from —0.1 to —0.3. On the other hand, 
when A goes larger, say A = 20 GeV, \J E ^ K is much larger than \/E — K and |ini(^-)p 
(|iin^-)P) is getting larger (smaller) as cos^^f changing from —0.1 to —0.3. It thus leads to 
the intersections for those curves with negative cos Qf values. 

Fig.s |^(a) and |^(b) show the mass ratio dependence of the tau polarizations. As 
shown in Fig. ^(a), -P(t'^) is a slow varying function of with fixed A value (A = 20 GeV 
in the plot). However, the ratio P{t^) / P{t^) varies significantly for 0.1 ^ cos6f ^ 0.3 with 
fixed A value as shown in 0(b). This plot may help in determining the model parameters 
such as and cos^^f in the Bino-like lightest neutralino scenarios. For example, a ratio 
P{t+)/P{t-) = -1.5 will lead to cos6'^ ~ 0.3 and > 2.5 in the model with A = 20 GeV. 

In Fig.s ^(a-c) we plot the tau polarization functions and their ratio as functions of 
cos^f for = 2 and five choices A = 4, 8, 12, 16, and 20 GeV. Fig. §(a) shows that the 
polarization -P(t+) increases as cos6'f increases for the cos6'f range under consideration. 
Similar behaviors occur when varying the A value in the model. Usually a larger A value 
will imply a lower tau polarization for (e^ —>■ e"^r+ff ) and a higher tau polarization for 
(e^ —>■ e^T^f^). Several curves for the tau polarization ratio P{t^)/P{t^) are plotted 
against the stau mixing parameter cos6'f with different A values in Fig. §(b). Varying the 

parameter could dramatically change the cos^f dependence of P{r^) / P{t^). We plot 
P{t^) / P{t^) for five choices = 1.0, 1.5, 2.0, 2.5 and 3.0 in Fig. §(c). As shown in the plot, 
the polarization ratio P{t^)/ P{t^) depends on r^j sensitively when 0.3 ^ | cos^^^ | ^ 0.2. 

We can also study the neutralino mixing effects on the tau-lepton polarizations. Fig. |^ 
shows typical plots in both the bino-like lightest neutralino scenario and in the Higgsino- 
mixed neutralino region with the model parameters mg^ = 200, \^\ = 300, mi = 400, 1712 = 
800 GeV and tan/3 = 10. By comparing Fig. I to |( a), we find that the tau polarization 
behavior in the Higgsino-mixed neutralino scenarios may look similar to that in the bino-like 
lightest neutralino scenarios with a larger | cos 6f\ value. For instance, the tau polarization in 
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the Higgsino- mixed neutralino region with cos 6f = —0.1 would look like the tau polarization 
in the bino-dominant lightest neutralino scenario with cos 6*^ ~ —0.3. 

In Fig. ^ we plot the normalized distributions for the final state electrons and tau leptons 
as functions of their transverse momenta pt in the selectron rest frame. Here we choose the 
model parameters mg^ = 200 GeV , = 1.2, cos^^ = —0.15. We plot three different 
values of A = 12, 16, 20 GeV respectively. As seen from the figure, the peak values occur 
at about pr = 7.5,5.8,4.3 GeV ( 6.0, 4.9, 3.7 GeV) for the final state tau lepton (electron) 
distributions as we vary the A from 20, 16, to 12 GeV accordingly. It should be noticed 
that when the mass difference A is too small, the final state tau may be too soft to detect 
and the determination of the tau polarization would be a challenge. However, if the mass 
splitting is large enough, the polarization of the final state tau lepton can be measured. 

According to the numerical results the knowledge of the tau polarization in the three- 
body slepton decays can provide another observable and may help to reveal the stau mixing 
angle and/or the bino-slepton mass ratio m^/mg^ in the small mass difference A << mg^j 
limit. Experimentally, there may be a way for measuring the polarized tau-lepton pro- 
ductions in the selectron decay by observing its subsequent tau decays ti^u^. and 
r+ tt^Ut-. The tau polarizations can be obtained by measuring the spectrum of £'^±/£'t-± 
in the decays. 

IV. CONCLUSION 

In this paper, we investigated the tau polarizations for the three-body slepton decays 
under the stau NLSP scenarios in the small mass difference limit A << rrig^. The exact 
formulas for calculating the slepton decay rates at the tree level are provided in section two. 

The measurements of tau lepton polarizations in the slepton decays provide another 
window for studying the SUSY model parameters. In the right-handed f/j limit, the tau 
asymmetry ratio Pr is solely determined by the mass difference parameter A. Numerical 
studies show explicit dependence on the stau mixing angles and/or the mass ratio m^/mg^. 
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We also studied the distributions of the final state electrons and tau leptons as functions 
of their transverse momenta pt in the selectron rest frame. The energy spectrum of the 
final state leptons are mainly determined by the mass difference ratio A = — ruf-^ . In 
our chosen parameters mg^ = 200 GeV, = 1.2, cos^f = —0.15 as we vary A from 12, 16 
to 20 GeV, the peaks of the distributions range from pt — 4.3, 5.8 to 7.5 GeV respectively. 
However if the mass difference is too small, the final state tau lepton will to too soft and 
the measurement of the polarization would be very difficult. But if A is large enough such 
that we can measure its polarization, the knowledge of tau-lepton polarization may help in 
determining the parameter space of a typical GMSB model. 
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FIGURES 

FIG. 1. The selection boost direction {(f), ip) in the lab frame where the electron moving direction 
and the tau moving direction spans the XZ plane. Here 6 denotes the angle between the outgoing 
electron and the outgoing tau lepton. 

FIG. 2. Tau polarization as the function of A in the three-body decays of cr with = 2 and 
cos^f =-0.3, -0.2, -0.1, 0.1, 0.2, 0.3. (a) P^(et ^ e+r+ff ) as function of A. (b) The P(r+)/P(r-) 
ratio as function of A. 

FIG. 3. Tau polarization as the function of the mass ratio = m^/rrie^ with seven choices 
cos6'f = -0.3, -0.2, -0.1, 0.1, 0.2, 0.3 and A = 20 GeV. (a) Pr(et e+r+ff) as function of r^. 
(b) The P{t^)IP{t~) ratio as function of r^. 

FIG. 4. (a) Tau polarization as function of cos^f for = 2 and five choices A = 4, 8, 12, 
16 and 20 GeV for the selectron decay e+r+ff. (b) The ratio P{t^)/P{t~) as function of 

cos6'f for = 2 and five choices A = 4, 8, 12, 16 and 20 GeV. (c)The ratio P(r+)/P(r-) for A 
= 20 GeV and five choices = 1.0, 1.5, 2.0, 2.5, 3.0. 

FIG. 5. The dependence of tau polarization on the mass difference A = mg^— mf^ with different 
cos Of values. Curves are plotted both in the bino-like neutralino scenario and in the higgisino-like 
neutralino scenario for the decay process e^ — > e+r+ff. Here we choose the model parameters 
mg^ = 200, \n\ = 300, mi = 400, = 800 GeV and tan/3 = 10. 

FIG. 6. Lepton px distributions in the rest frame of the e^ decaying to e+r+f ~. Normalized 
Pt distributions for final e (dot-dash lines) and r (solid lines) for the model parameters with 
m-e^ = 200 GeV, cosfilf = -0.15, = 1.2, and A = 12, 16, 20 GeV respectively. 
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